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ABSTRACT 

Infrared magnitude-redshift relations (K-z relations) for the SCR and 6C samples 
of radio galaxies are presented for a wide range of plausible cosmological models, 
' including those with non-zero cosmological constant Q\. Variations in the galaxy 

formation redshift, metallicity and star formation history are also considered. The 

■ results of the modelling are displayed in terms of magnitude differences between the 
^ ' models and no-evolution tracks, illustrating the amount of K-band evolution necessary 

I to account for the observational data. 

■ Given a number of plausible assumptions, the results of these analyses suggest 
' that: (i) cosmologies which predict Tq x Hq > I (where Tq denotes the current age of 

the universe) can be excluded; (ii) the star formation redshift should lie in the redshift 
interval 5 < Zf < 20, values towards the lower end of the range being preferred in 
cosmologies with larger values of Tq x Hq; (iii) the Einstein-de Sitter model (Oq = 
, 1, fl\ = 0) provides a reasonable fit to the data; (iv) models with finite values of 

' can provide good agreement with the observations only if appropriate adjustments 

of other parameters such as the galaxy metallicities and star-formation histories are 
Q . made. Without such modifications, even after accounting for stellar evolution, the high 

i—{ ' redshift radio galaxies are more luminous (that is, more massive) than those nearby 

I in models with finite Qa, including the favoured model with — 0.3, = 0.7. 

d • For cosmological models with larger values of Tq x Hq, the conclusions are the same 

regardless of whether any adjustments are made or not. The implications of these 
results for cosmology and models of galaxy formation arc discussed. 
■ 

■ Key words: galaxies - active galaxies - radio galaxies - galaxy evolution - cosmo- 
' logical parameters. 
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1 INTRODUCTION 

The models used to account for the infrared magnitude- 
redshift relation of luminous radio galaxies are sensitive to 
the assumptions made concerning the evolution of their stel- 
lar content, as well as the choice of cosmological model. 
Galaxies at a particular redshift are observed at different 
cosmic epochs in different cosmological models and so their 
predicted spectra at a given redshift can vary quite dramat- 
ically. Accurate modelling of the evolution of galaxy spectra 
with cosmic time is required if the infrared K magnitude- 
redshift relation is to be used to constrain the range of ac- 
ceptable world models. In this paper the spectral evolution 
of galaxies is modelled for different plausible cosmological 
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models using the most recent spectral synthesis codes of 
Bruzual and Chariot (2001), GISSEL2000, and these results 
compared with the observational data. Even for the sim- 
plest cases, numerous cosmological and evolutionary effects 
need to be taken into account and it is often not intuitively 
obvious how changes in the assumptions made change the 
predicted relations - one of the objectives of this paper is to 
provide these insights. Although our prime concern is the K- 
z relation for radio galaxies, the considerations of this paper 
are relevant to all such cosmological studies. 

The radio galaxies associated with the most powerful 
extragalactic radio sources seem to form a remarkably uni- 
form class of galaxy out to large redshifts. Following the 
early work of Lilly and Longair (1984), Best et al. (1998) 
carried out a much more detailed analysis of the photomet- 
ric properties of a complete sample of SCR radio galaxies, 
the brightest radio galaxies in the northern sky selected at 
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178 MHz (Laing et al. 1983). They modelled the K-z relation 
for world models with f2o = and flo = 1, in both cases as- 
suming Qa = 0. These models included passive evolution of 
the stellar populations of the galaxies with formation red- 
shifts of 3, 5 and 20. The passive evolution models in the 
Clo = 1 universe provided a reasonable fit to the observa- 
tions, whereas those with f2o = gave a poor fit. The mod- 
els with low redshifts of formation, Zf ~ 3, resulted in more 
luminous galaxies than those observed at z ^ 1.2, suggest- 
ing that the stellar populations of the galaxies must have 
formed at z > 3. 

Bales and Rawlings (1996) and Bales et al. (1997) de- 
termined the K-z relation for a complete sample of 6C radio 
galaxies which were selected at 151 MHz and are about six 
times fainter in fiux density than the 3CR sample (Bales 
1985). The 6C host galaxies have similar K-luminosities to 
those of the SCR galaxies at z < 0.8, but are about 0.6 
magnitudes fainter at larger redshifts. Nevertheless, the high 
redshift 6C galaxies are significantly brighter than expected 
in world models in which no evolution of the stellar popula- 
tions of the galaxies is assumed. 

The objectives of this paper are to review the current 
status of the K-z relation for 3CR and 6C luminous radio 
galaxies (Section 2) and then to carry out a more extensive 
analysis of the resulting K-z relation for a wider range of 
plausible world models than those considered by Best et al. 
(1998) and Bales et al. (1997), including models with a non- 
zero cosmological constant. The predicted evolution of the 
stellar populations of the galaxies is derived from the most 
recent spectral synthesis codes of Bruzual & Chariot (1993, 
2001), GISSEL2000; a wide range of stellar synthesis models 
is considered, as well as the effects of changes in metallicity. 
The details of the different cosmological models, the spec- 
tral synthesis codes and procedures used to predict the K-z 
relations are described in Section 3. The results of these 
calculations are described in Section 4, indicating which pa- 
rameters have the greatest impact upon the form of the pre- 
dicted K-z relation. Also included in this Section are the 
results of varying the star formation history and metallicity 
in plausible ways. Most massive galaxies at small redshifts 
have high metallicities and colour-magnitude relations show 
that more massive galaxies are redder in colour, suggesting 
a mass-metallicity correlation (see, for example, Arimoto & 
Kodama, 1996). On the other hand, at high redshifts, we 
might expect a lower metal content, as suggested by the 
analyses of Pei et al. (1999) and Pettini et al. (1997). In 
Section 5 the results of these analyses are discussed and our 
conclusions are summarised in Section 6. Throughout the 
paper. Ho = 50 km s~^ Mpc~^ is adopted, unless otherwise 
stated. 



2 THE K-Z RELATION FOR LUMINOUS 
RADIO GALAXIES 

2.1 The Revised K-z Diagram 

The observational data used in this analysis consist of in- 
frared K magnitudes for 72 3CR galaxies in the redshift 
interval < z < 1.8 (Best et al. 1998), and 57 6C galax- 
ies in the range < z < 3.4 (Bales et al. 1996, 1997). The 
3CR galaxies comprise a complete sample of radio sources 



with fiux densities Sns ^ 10.9 Jy at 178 MHz, with a mean 
flux density of about 15 Jy; the 6C sample is complete in 
the flux density interval 2.2 < S151 < 4.4 Jy at 151 MHz, 
with an equivalent mean fiux density at 178 MHz of about 
2.6 Jy. The 3CR K magnitudes were obtained for z < 0.6 
from observations with the IRCAM3 array on the UK In- 
frared Telescope (Best et al. 1997), and for z < 0.6 from the 
photometry of Lilly & Longair (1984). The 6C data points 
were mostly derived from IRC AM observations, but also from 
colour corrected K' observations with the redeye camera 
of the Canada- France-Hawaii Telescope (Bales et al. 1997). 

Various factors had to be included to create a self- 
consistent set of K magnitudes. Firstly, aperture corrections 
were made, reducing the magnitudes to a standard circular 
aperture of metric diameter 63.9 kpc. These aperture cor- 
rections are dependent on the cosmological model, due to 
differences in the angular diameter versus redshift relation; 
however, the changes from the standard Oq = 1, Oa ~ 
model are generally small, ^0.1 mag, and so the data were 
converted to 63.9 kpc in the Slo ~ 1, = model, and the 
model-to-model variations were included in the predicted 
model tracks. 

The 63.9 kpc aperture was chosen since this had been 
adopted for the 6C observations wherever possible. In cases 
where there was a nearby companion galaxy, or other prob- 
lems prevented the use of the full aperture, magnitudes for 
the 6C galaxies were determined for smaller apertures and 
aperture corrections applied by Bales et al. (1997); for galax- 
ies at redshifts z < 0.6 they assumed that the radio galaxies 
follow an elliptical galaxy growth curve (Sandage 1972), and 
at redshifts z > 0.6 they assumed that the growth curve for 
the luminosity L of the galaxy within radius r could be rep- 
resented by L(< r) oc r"'^^. The K-magnitudes of the 3CR 
galaxies were originally determined within a 9 arcsec diam- 
eter aperture. At large redshifts, this is almost equivalent to 
63.9 kpc, and so the aperture corrections were small. Larger 
aperture corrections were necessary at low redshifts since 
the 9 arcsec aperture included only the central regions of the 
galaxy (Lilly & Longair 1984). A de Vaucouleurs profile for 
giant ellipticals was used to correct these magnitudes to the 
standard diameter, since the analysis of Best et al. (1998) 
showed that this is a good representation of the variation 
of the K surface brightness with radius. It is also consistent 
with the radial optical light profile found for low redshift 
radio galaxies by Lilly, McLean and Longair (1984). 

The resulting K-z relation for 3CR and 6C galaxies, rep- 
resented by open boxes and crosses respectively, is shown 
in Fig. ^. New versions of the predicted K-z relations are 
shown for an f2o ~ 1, ^a = world model. These involve 
the use of the most recent set of spectral synthesis codes 
(gissel2000), for which a Scalo Initial Mass Function was 
adopted, with upper and lower mass cut-offs of IOOM0 and 
O.IM0 respectively; solar metallicity was assumed. In ad- 
dition, improved K-band transmission data for the K-band 
filter of UKIRT were used. Details of these procedures are 
discussed in Section 3.2. The differences as compared with 
the analysis carried out by Best et al. (1998) are very small, 
and so their conclusions for the f2o = and fio = 1 models 
are unchanged. In particular, the flo — 1, Qa = model 
with passive evolution of the stellar content of the galaxies 
provides a satisfactory fit to the observations and will be 
used as a fiducial model in the new computations. 
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Figure 1. The revised K-z diagram for powerful radio galaxies for the world model with f2o = li^^A = 0. The open squares represent 
SCR galaxies and the crosses 6C galaxies. The solid curve corresponds to the predicted K-z relation for a non-evolving stellar population. 
The three other relations are for passively evolving galaxies in which the stellar population formed at a constant rate for an initial period 
of 1 Gyr. The three formation redshifts were: 2f = 3 - dot-dashed line; = 5 - dashed line, and Zf = 20 — dotted line. The solid triangles 
show mean values of the K-magnitudes for both the SCR and 6C data averaged in intervals of 0.1 To, where To is the age of the world 
model, To = 2/3_ffo for this case. The bars show the standard deviation of the distribution of magnitudes in each of these bins. 



2.2 The SCR, 6C and MRC Radio Galaxies 

One of the intriguing features of the comparison of the SCR 
and 6C data is that, although they are in excellent agree- 
ment at redshifts z < 0.6, there is a progressive shift in 
their mean magnitudes at greater redshifts, the mean K- 
magnitude of the 6C galaxies at 2; > 1 being about 0.6 
magnitudes fainter on average than the SCR galaxies (Eales 
et al. 1997). This can be seen in Fig. where the boxes, 
representing SCR galaxies, lie towards the low limit of the 
distribution of the crosses, 6C galaxies, at 2 > 1. 

Given the quite large dispersion about the mean K-z re- 
lation, Bayesian statistics have been used to investigate the 
significance of this magnitude difference and how it evolves 
with redshift without making restrictive assumptions about 
the underlying distribution (see Figure |^ . The results show 
clearly the mean difference between the two data sets of 
about 0.6 magnitudes at high redshift /early times and its 
decrease to zero at later times. The mean value of 5y(x) 
increases slightly at the lowest redshifts to about 0.15 mag- 
nitudes, but this increase is not significant. At high redshifts 
the difference in K-magnitudes between SCR and 6C galax- 
ies is significant at the 95% confidence limit. 

This result for the 6C galaxies is mirrored by that of the 
still fainter 7C radio galaxies (Willott, priv comm). However, 
McCarthy (1999) and van Breugel (1998) have studied the 
K-z relation for radio galaxies in the MRC IJy sample, which 
contains sources with flux densities 5408 > Uy at 408MHz 



(= Si78 ~ 5.8Jy) from the Molonglo reference catalogue 
(McCarthy et al 1996), and they find that the mean K- 
band magnitudes of the MRC radio galaxies to be only 
to 0.2 magnitudes fainter than the SCR, a much smaller 
difference than for the 6C sources. The ~ 25% spectroscopic 
incompleteness of this sample may introduce a bias in that 
it is more difficult to obtain redshifts for intrinsically fainter 
galaxies, but McCarthy argues that this bias is unlikely to be 
sufficiently great to alter the mean magnitude of the MRC 
galaxies by 0.6 magitudes. 

The important aspect of this discussion is whether or 
not there is a correlation between radio luminosity and K- 
absolute magnitude, as would be the case if it is accepted 
that there is a real difference between the K-z diagrams for 
the SCR and 6C samples. For the purpose of the present 
analysis, we have assumed that the difference between the 
6C and SCR samples, being smaller than the scatter within 
either sample, can be ignored. Therefore, we have taken 
mean values and dispersions for the combined sample. These 
are indicated by the solid triangles and associated error bars 
in Fig. I 

A literal interpretation of the K-z diagram shown in 
Fig. ^ suggests that to a good approximation the luminous 
radio galaxies belong to a single population of passively 
evolving galaxies of the same stellar mass. The galaxies also 
have the same characteristic radii of about lOkpc over the 
entire range < z < 1.5 (Best et al 1998, McLure & Dun- 
lop 2000), which supports this conclusion. As discussed by 
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Figure 2. The inferred K-band magnitude difference (dy(x)) be- 
tween the 3CR and 6C galaxy populations as a function of cosmic 
time, expressed as a fraction of Tq, the present age of the universe, 
determined using a Bayesian statistical approach. In the Bayesian 
approach, two flexible distributions for the reduced magnitude 
y{x) (the difference between the observed magnitudes K and the 
expectation Kq of the no evolution model) are assumed to under- 
lie the two data sets, for the range of redshifts covered by both 
samples. y{x) is represented by a sum of Legendre polynomials 
such that y(x) = Tii^Wf^^f^lx). Similarly, the scatter in the ob- 
servational data (j{x) is represented by log^g (t{x) = llbh^i^{x). 
The differences between the two data sets and their scatter are 
defined by the quantities Sy{x) = S(5«)(,#^(x), and <51og]^Q a-{x) = 
T,Sbh^h{x) respectively. Gaussian priors were assigned to the pa- 
rameters tOfj, feft, Swf^ and 56^; the widths of these four Gaussians 
were hyperparameters inferred from the data using the BUGS soft- 
ware (Bayesian inference Using Gibbs Sampling; Thomas et al. 
1992). This method is similar to that used by MacKay (1995). The 
best estimate of the difference in magnitude between the samples 
with cosmic time is plotted as the central solid line. The 95% 
confidence limits are represented by the shaded region, delimited 
by two solid lines. The cosmological model used has Qq = 0.3 
and fl\ = 0.7. 



Best et al, this is a somewhat surprising result, since the 
preferred hierarchical clustering picture of galaxy formation 
predicts that the masses and radii of galaxies would grow 
significantly over the redshift range considered. Evidence 
suggests that galaxy environments also differ with redshift; 
high redshift SCR galaxies seem to belong to proto-cluster 
environments, whereas low redshift galaxies are more com- 
monly found in small groups or otherwise isolated environ- 
ments. Intuitively, it seems difficult to reconcile these points 
with closed-box passive evolution of the stellar populations 
of the SCR galaxies. Best et al. recognised this problem and 
argued that the luminous radio galaxies must be observed at 
a similar stage in their evolution, specifically that at which 
they had attained a certain stellar mass. Provided that stars 
brought together by on-going mergers were formed at the 
same early epoch, the old stellar populations of galaxies 
formed by hierarchical clustering would appear to be pas- 
sively evolving. The small scatter in the colour-magnitude 
relations of elliptical galaxies observed out to redshifts 2 ~ 1 
(Stanford et al. 1998, Bower et al. 1998) justifies this as- 
sumption. 

A further objective of our analysis is to determine 
whether it remains true that the luminous galaxies present 



on the K-z relation have the same stellar mass, and how dif- 
ferent assumptions about the underlying cosmology and the 
evolution of the stellar populations change this conclusion. 



3 MODELLING PROCEDURES 
3.1 Cosmological Models 

Because of the availability of the GISSEL codes of Bruzual 
and Chariot, the opportunity has been taken to consider 
a wider range of evolutionary histories of the star forma- 
tion rate of radio galaxies and to couple this with an ex- 
panded range of cosmological models, in particular, those 
with non-zero values of the cosmological constant (A) . Fried- 
man world models with zero cosmological constant and 
Qo — 0, 0.1, O.S and 1.0 have been considered, as have a range 
of models with non-zero Oa. The relation between cosmic 
time and redshift is: 



Jo Ho J, (l + 2)[r!o(l + 2)3 + OA]i/2- 

For the world models with flat geometries, Qo + Qa ~ 1, 
corresponding to those favoured by recent analyses of the 
power spectrum of fluctuations in the Cosmic Microwave 
Background Radiation as observed in the Boomerang and 
Maxima experiments (e.g. Jaffe et al. 2000), there is a para- 
metric solution: 

2 , /l + cose\ , fQoV^^,.. ^3/ 

(2) 

Numerical integrations have been used for other models, and 
for determining the variation of co-moving radial distance 
coordinate with redshift. The models with non-zero SIa con- 
sidered in this analysis are listed in Table 1. Two of these 
have flat geometries and the other two models have open 
geometries with low values of both f2o and Q,h. This range 
of models with flnite values of Qa spans those which would 
be consistent with the results found in the Supernova Cos- 
mology project (Perlmutter et al. 1999). 

3.2 Galaxy Spectral Synthesis 

The GISSEL codes of Bruzual and Chariot provide a means of 
modelling galaxy spectra through isochrone synthesis. The 
models used in the present analysis all include the passive 
evolution of the stellar populations of the galaxies. The stan- 
dard model is taken to have solar metallicity, Z = Zq, and 
the stellar population of the galaxy was formed in an initial 
uniform starburst with a duration of 1.0 Gyr. The GISSEL 
code enabled the spectral evolution of this model to be de- 
termined for ages up to 20 Gyr and this is shown in Fig. ^(a) 
at a number of ages which are equally spaced logarithmically 
in cosmic time. Generally, only the late stages of evolution of 
the starburst (t 3> 10^ years) are of interest in this analysis. 
For one of the models, Q.o =0.1 and Oa = 0.9, the present 
age of the Universe, To, is greater than 20 Gyr and evolu- 
tion in this model was terminated at a low redshift limit of 
z = O.S. 

It is instructive to note that, in Fig. S(a), at optical 
and infrared wavelengths the spectra are roughly equally 
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Table 1. Parameters of the different cosmologies investigated in 
this paper. It is assumed that Hq = 50 km s~^. Mpc~^ 





Ha 


Age at ij = 1 
(To X Ho) To (Gyr) 


1.0 


0.0 


0.67 


13.0 


0.3 


0.0 


0.81 


15.8 


0.1 


0.0 


0.90 


17.6 


0.0 


0.0 


1.00 


19.6 


0.3 


0.7 


0.96 


18.9 


0.3 


0.3 


0.86 


16.9 


0.1 


0.9 


1.28 


25.0 


0.1 


0.3 


0.98 


19.1 



spaced logarithmically in luminosity at ages greater than 
10^ years, indicating that, at a given wavelength, the spec- 
tral energy decreases roughly as a power-law with cosmic 
time. As shown by Gunn (1978), Longair (1998) and others, 
this power-law behaviour can be understood in terms of the 
rate at which stars evolve from the main sequence onto the 
red giant branch, if the initial mass function is assumed to 
be of power-law form. Such simple arguments predict that 
L cc for the f^o = 1, ^^a ~ model. This is not so dif- 
ferent from the evolution seen in Fig. 3(a), but there are 
important variations in the power-law index as a function 
of wavelength. At 0.55 /im, the evolution can be roughly 
described by L0.5 oc t^^'^, whereas at 2.2 fim the exponent 
decreases to about —0.5. Therefore, the simple approxima- 
tion is not adequate for the present analysis. For essentially 
all our computations, only evolution at epochs greater than 
lO'' years is required. It can be seen however that at ear- 
lier times the K-band luminosity increases with age. This is 
associated with the evolution of stars onto the asymptotic 
giant branch. The stellar evolution models of Maeder and 
Megnet (1989) and Girardi et al. (2000) show that, for stars 
with M > 2M0, as the stellar mass decreases, the luminos- 
ity of the upper AGB stars increases and their temperature 
decreases. Thus, the simple monotonically decreasing spec- 
tral luminosity is only applicable for times greater than 10® 
years when the turn-off mass drops below 2 M0 and the AGB 
contribution becomes less important than evolution onto the 
red giant branch. 

In previous analyses, it was assumed that the metallic- 
ity of the galaxies corresponded to solar values. The evidence 
of the Madau plot in its various incarnations (see, for exam- 
ple, Madau et al. 1998) and the analyses of the metallicities 
of the absorption line systems in distant quasars have sug- 
gested that the chemical abundances of the elements have 
evolved very significantly over the redshift interval < z < 4 
(see, for example, Pettini et al. 1998). To understand the 
effects of metallicity upon the predicted K-z relation, mod- 
els have been created for four different metallicities: 0.2^0, 
0.4^0, Zq, and 2.5^©. 

Different star formation histories have also been consid- 
ered. These include a continuous exponentially decreasing 
star formation rate, with an e-folding time of 1 Gyr, and 
a simple delta-function starburst. Combined models could 
then be constructed involving multiple star formation events 
occurring at different cosmic times. For example, if the on- 
set of radio source and quasar activity were contemporane- 



ous with (or stimulated) an active starburst, the additional 
starburst component might affect the K-magnitude of the 
galaxy. To illustrate how such events could impact the K-z 
relation, model starbursts can be added to the underlying 
passive evolution of the majority old stellar population. Two 
model star bursts were considered. In the first, the starburst 
took place instantaneously, and was observed at times 10^, 
10^ and 10* years after that event; this spectral evolution is 
shown in Fig. ^(b). The intriguing aspect of these spectra 
is that the starburst is most luminous in the infrared wave- 
band after about 10^ years. After about this time the most 
massive stars complete their evolution on the main sequence 
and join the asymptotic giant branch. It is interesting that 
the time scale of 10^ years is of the same order as the typical 
age of the radio source events associated with the SCR and 
6C sources and so such a star burst could contribute to the 
K-luminosity of the galaxy. In the second example, the case 
of a constant star burst lasting 10* years is considered, but 
observed at 10^, lO'^ and 10* years after the start of that 
event; these spectra are shown in Fig. ^(c). In this case, the 
star burst is expected to have a much smaller effect upon 
the K-luminosity of the galaxy. It can be seen that, after 
lO'' years, the starburst has the same spectrum as the first 
model, but with only one hundredth of the luminosity since 
only that fraction of the stars have formed by that time. 
On the other hand, after 10* years, the starburst is more 
luminous than in the first model since there is a significant 
contribution of recently formed stars to the K-luminosity. In 
either case, it can be seen that starbursts observed between 
10'' and 10* years after the onset of the starburst can have 
a significant impact upon the K-luminosity of the galaxy. 
Models have been considered in which the mass fractions 
of newly formed stars range between 0.5% and 5.0% of the 
total stellar mass of the galaxy. 

The resulting spectral energy distributions (SEDs) were 
then redshifted and convolved with the transmission func- 
tion of the standard K-filter. These were then combined with 
the different world models to create K-z diagrams. A cata- 
logue of the range of models considered is given in Tables 3 
and 4. 



4 RESULTS 

4.1 Changing the Cosmological Model 

It is simplest to illustrate the effects of changing the parame- 
ters of the cosmological models by subtracting the predicted 
K-z track for the chosen cosmology from the corresponding 
predictions of the flo = 1.0, Qa = model: AK is defined 
as K(Slo = 1) — K(model) so that negative values of AK 
mean that a galaxy of a given absolute luminosity would be 
observed to have a fainter apparent magnitude. The results 
for the cosmologies listed in Table 1 are shown in Fig. ^ and 
tabulated in Table 2. 

It is instructive to distinguish those contributions to 
AK due to the choice of cosmological model and those due 
to the evolving spectrum of the galaxy. The latter can be 
thought of as the K-correction including the effects of the 
evolution of the stellar population. Following the notation 
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Figure 3. (a) Synthesised Spectral Energy Distributions for the 
standard passive evolution model, consisting of an initial uniform 
starburst of duration 10* years with solar metallicity, Z = Zq. 
The predicted spectral energy distributions of the galaxy is dis- 
played at eight different ages which are equally spaced logarith- 
mically in cosmic time. The vertical lines show the regions of the 
spectrum corresponding to the rest wavelength of the observed 
K-band radiation at redshifts z = 3, 2, 1 and 0. (b) SEDs for 
the additional starbursts used in our alternative star formation 
histories, as described in section 3.2. These are instantaneous star- 
bursts, observed at ages of 10®, 10^ and 10* years. 



Figure 3. — Continued, (c) The spectral evolution of a uni- 
form 10* year starburst observed at ages 10^, 10^ and 10* years, 
the spectrum at the latest time being the brightest over the full 
range of wavelength plotted. In all three figures, the spectra are 
normalised to the luminosity per unit solar mass of stars. 



of Longair (1998), the flux density of the galaxy is given by 
L[vq{1 + z),U{z)] 



Si{vo, z) = 



(3) 



where i labels the choice of cosmological model and Di{z) 
is a distance measure which is related to the standard lumi- 
nosity distance Dl by Dl — Di{z){l + z). The luminosities 
of the modelled galaxies in the different cosmologies are nor- 
malised by varying their stellar mass, such that at zero red- 
shift they have K-band luminosities -Lk(z=o) matching the 
observed K-z relation. Therefore if L'[v,t] is the luminosity 
of the galaxy per solar mass, 

L[vo{l + z)Mz)]^L'[uo{l + z)Mz)] X ,,t.'^^T°i,i (4) 



The ratio of the flux densities in the ith model to that in the 
fiducial model, for the same observed flux density at small 
redshifts is given by 

Dljz) 
(5) 

Thus, the magnitude difference associated with the choice 
of cosmological model is 



Si{uo,z) ^ L'[uo{l + z),U{z)] 



L'Wo,to{0)] 
L'[uo{l + z),to{z)] 



AKn = 51og[A(^)/A)(^)l, 



(6) 



while that associated with observing different parts of the 
evolving spectral energy distribution at different cosmic 
times is 



AKov = 2.5 log 



L'[uo{l + z),U{z)] 



L'[uo,to{0)] 



L'[uo{l + z),to{z)] 



(7) 



The total change in magnitude is given by AK = AKn + 
AKev Table 2 shows these differences for the models listed 
in Table 1. 

Fig. 4 shows that galaxies are expected to be fainter 
at large redshifts in all the models considered as compared 
with the model with Slo = 1, f^A = 0; this means that for 
the same observed K-magnitude, the galaxies would have 
to be intrinsically more luminous. There are several effects 
which contribute to this decrease. The most important of 
these is the fact that, for a given redshift, galaxies are ob- 
served at greater radial co-moving coordinate distances as 
compared with the model with f2o = 1, JIa = 0. In addition, 
the cosmic times at which the spectral energy distributions 
are observed are significantly different. What is important, 
however, are the relative epochs at which the galaxies are 
observed. It turns out that, in the passively evolving mod- 
els, the relative changes in observed luminosity (AKcv) are 
significantly less than the effects of the cosmological model 
(AKn). This can be understood from Table 2 which shows 
the ages and distance coordinates at which galaxies at differ- 
ent redshifts are observed in the different cosmologies listed 
in Table 1 assuming Ho = 50 km s~^ Mpc"'^; included in 
this Table are the contributions of AKn, the total change in 
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Table 2. The cosmic time and distance coordinate D, in units of c/Hq at which galaxies are observed in the cosmological models listed 
in Table 2 at redshifts 0, 1, 2 and 3. It is assumed that Hq = 50 km s~^ Mpc~^. The third line of each entry shows the contribution of 
the choice of cosmological model AKq to AK. The fourth line shows the exact contribution of AKcv to AK. The fifth line shows the 
total value of AK which is derived from integration over the band pass of the K-filter. 



Ho = 1 
Ha = 



no = 



0.3 
= 



f^o = 

f^A 



0.1 
= 



no ■■ 



: 0.3 

= 0.7 



Qa 



■ 0.3 
= 0.3 



Ho ■■ 

flA 



: 0.1 

= 0.9 



Qa 



0.1 

: 0.3 



13.0 Gyr 
D =0.00 
AKn=0.00 

AKev=0.00 

AK =0.00 



15.8 Gyr 
D = 0.00 
AKn= 0.00 
AKev= 0.00 
AK = 0.00 



17.6 Gyr 
D = 0.00 
AKn= 0.00 
AKcv= 0.00 
AK = 0.00 



19.6 Gyr 
D = 0.00 
AKn= 0.00 
AKov= 0.00 
AK = 0.00 



18.9 Gyr 
D = 0.00 
AKn= 0.00 
AKev= 0.00 
AK = 0.00 



16.9 Gyr 
D = 0.00 
AKn= 0.00 
AKcv= 0.00 
AK = 0.00 



25.0 Gyr 
D = 0.00 
AKn= 0.00 
AKev= 0.00 
AK = 0.00 



19.1 Gyr 
D = 0.00 
AKq= 0.00 
AKev= 0.00 
AK = 0.00 



4.61 Gyr 
D =0.59 
AKn=0.00 
AKcv=0.00 
AK =0.00 



6.54 Gyr 
D = 0.69 
AKn=-0.34 

AKev = -0.13 

AK =-0.47 



7.97 Gyr 
D = 1.11 
AKn=-0.47 
AKcv= 0.00 
AK =-0.47 



9.80 Gyr 
D = 0.75 
AKn=-0.54 

AKev = -0.11 

AK =-0.65 



8.05 Gyr 
D = 0.77 
AKn=-0.60 
AKev=-0.02 
AK =-0.62 



7.07 Gyr 
D = 0.72 
AKn=-0.44 
AKcv=-0.04 
AK =-0.48 



12.73 Gyr 
D = 0.90 
AKn=-0.92 

AKev = -0.21 

AK =-1.13 



8.89 
D -- 
AKq = 

AKev = 

AK = 



Gyr 
= 0.77 
=-0.59 
=-0.07 
=-0.66 



2.51 Gyr 
D =0.85 
AKn=0.00 

AKev=0.00 

AK =0.00 



4.34 Gyr 
D = 1.10 
AKn=-0.59 
AKcv=-0.29 
AK =-0.88 



4.92 Gyr 
D = 1.24 
AKn=-0.84 
AKov=-0.13 
AK =-0.97 



6.53 Gyr 
D = 1.33 
AKn=-0.99 
AKcv=-0.39 
AK =-1.38 



4.52 Gyr 
D = 1.12 
AKn=-0.78 
AKcv=-0.10 
AK =-0.88 



4.07 Gyr 
D = 1.15 
AKn=-0.67 
AKov=-0.12 
AK =-0.79 



7.55 Gyr 
D = 1.54 
AKn=-1.30 
AKev=-0.40 
AK =-1.70 



5.45 Gyr 
D = 1.32 
AKq=-0.97 
AKov=-0.21 
AK =-1.18 



1.63 Gyr 
D =1.00 
AKn=0.00 

AKev=0.00 

AK =0.00 



2.57 Gyr 
D = 1.43 
AKq=-0.77 
AKev=-0.72 
AK =-1.49 



3.46 Gyr 
D = 1.68 
AKq=-1.13 
AKov=-0.42 
AK =-1.55 



4.90 Gyr 
D = 1.88 
AKn=-1.37 
AKcv=-0.85 
AK =-2.22 



2.96 Gyr 
D = 1.49 
AKn=-0.86 
AKev=-0.42 
AK =-1.28 



2.72 Gyr 
D = 1.46 
AKn=-0.82 
AKcv = -0.37 
AK =-1.19 



5.04 Gyr 
D = 1.98 
AKn=-1.48 
AKov=-0.81 
AK =-2.29 



3.80 Gyr 
D = 1.77 
AKf2=-1.24 
AKov=-0.39 
AK =-1.63 



K- magnitude AK and the difference between them AKcv, 
which represents the contribution of the changing spectral 
energy distribution of the galaxies with cosmic epoch. 

Table 2 demonstrates quantitatively the important 
point that, at all redshifts, the choice of cosmological model 
is always much more important than the evolution correc- 
tions within the context of passively evolving models for the 
stellar populations of the galaxies, so far as the differences 
from the fiducial model are concerned. 

It is interesting that it is possible to reproduce quite 
accurately the results of the exact computations by interpo- 
lating the luminosity changes at the appropriate redshifted 
wavelength in the template spectra shown in Fig. 3(a). It 
is found that this approach generally produces results accu- 
rate to about 0.1 mag, which is normally accurate enough 
to draw useful conclusions. The value of this approach is 
that it enables the investigator to use Fig. 3(a) to work out 
quite precisely the effects of the K-corrections for any other 
world model with any given value of Hubble's constant and 
for alternative assumptions relative to the template spectral 
evolution shown in Fig. 3(a). 

Poggianti (1997) provides tables of K-corections and 
evolutionary corrections for observations in a number of spe- 
cific filters, assuming an f7o = 0.45 cosmology. We have com- 
pared these results with those of our two most similar cos- 
mologies, i.e. fio = 1, ^^A = and flo = 0.3, f^A = 0. Despite 
a different set of spectral synthesis models, the evolutionary 
corrections for an elliptical galaxy formed with an e-folding 
time of 1 Gyr given by Poggianti are in good general agree- 
ment with our own, suggesting that the choice of stellar 



synthesis code does not have a large impact on the evolu- 
tionary corrections determined. The K-corrections given by 
Poggianti for their single cosmology in the K waveband also 
agree well with our own results. 

Considering first the range of models shown in Fig. 4, 
it is evident that only the passively evolving Slo = 1, f^A = 
model provides a reasonable fit to the data. The error bars 
shown in the figure represent the standard errors in the mean 
of the K-magnitude of the galaxies, within the given redshift 
intervals, and are less than the standard deviations shown 
in Fig. 1 by a factor of \/iV, where iV ~ 10 - 20. 

For illustrative purposes, a quantitative comparison be- 
tween four of the models plotted in Fig. ^ is given in Table 
3. These are the models with f2o = 1.0, Q,a ~ 0.0; flo — 
0.1, Qa = 0.0; ^0 = 0.1, fiA = 0.3; and D,o = 0.3, Oa = 0.7. 
In the first part of the table, the magnitude difference, AK, 
between the data points and the no-evolution tracks at red- 
shifts z — 0.5, 1.0 and 2.0 are listed, and in the second the 
difference in AK between the data points and the passive 
evolution models of Fig. ^ are shown. These data are also 
included as the first two entries in Table 4. The K-z rela- 
tions have been normalised to the same K-magnitude at low 
redshift, for both the passive evolution and the no-evolution 
tracks since the mean magnitudes of the galaxies at z < 0.3 
are consistent with a very well-defined linear K-z relation, 
as can be seen from Fig. 1. Positive values of AK indicate 
that the predicted magnitudes correspond to brighter galax- 
ies than are needed to account for the data. For example, the 
Qo = 1.0, f^A = model results in slightly brighter tracks 
than the mean of the data points at z = 2.0. 
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0.1 

Redshift 



Figure 4. Changes in the predicted K-z relations for different 
cosmological models. The predictions are presented as differences 
relative to that of the Hq = 1, S^a = model. AK is defined 
as K(Oo = 1) ~ K(model), so that negative values for AK mean 
that a galaxy of given absolute luminosity would be observed 
to have a fainter apparent magnitude. All the models include 
evolution of the stellar populations of the galaxies according to 
the standard passive evolution model. Models with Ha = and 
Qo = 0.0, 0.1, 0.3 and 1.0 are shown as dotted lines, the last model 
having AK = by definition; AK becomes more negative with 
decreasing Qq. Models with Qq = 0.1, Qa = 0.9 and Oo = 0.1, 
= 0.3 are shown by dashed lines, the latter model having the 
more negative value of AK. Models with Qq = 0.3, Qa = 0.3 
and Qq = 0.3, Ha = 0.7 are shown by solid lines, the latter 
model having the more negative value of AK. The triangles and 
associated error bars represent the observed values from Fig. hi. 



The importance of these results is that, in general, ex- 
cept for the fiducial model, Qo = 1,Qa = 0, the galaxies 
are predicted to be fainter at large redshifts than they are 
observed to be, even when passive evolution is taken into 
account. This is a somewhat unexpected result. First of all, 
according to the hierarchical clustering model, galaxies at 
large redshifts should be less luminous than those nearby 
since they contain fewer stars. The present results go in en- 
tirely the wrong direction. For the most extreme To x Ho 
models the galaxies would need to contain about ten times 
more stars to account for these differences. Second, it might 
be argued that this result would be found if there were a 
correlation between radio luminosity and the K-luminosity 
of the galaxy. In this case, the galaxies at large redshift 
would be expected to be larger than those at small red- 
shifts. The analysis of McLure and Dunlop (2000) showed, 
however, that the half-width diameters of the luminosity 
profiles of the radio galaxies are invariant as a function of 
redshift, suggesting that the galaxies have attained a simi- 
lar evolutionary state so far as their total stellar populations 
are concerned. 

Because of these somewhat surprising conclusions, the 
impact of changing the assumptions made in the fiducial 
model in plausible ways has been analysed in more detail. 




1 1 1 1 1 — 

0.5 



Cosmic Time as a fraction of To 




1 1 1 1 1 1 1 1 1 1 1 

0.5 1 

Cosmic Time as a fraction of To 
Figure 5. The values of AK relative to the no-evolution model as 
a function of cosmic time, for galaxies with different assumed val- 
ues of the redshifts, Zf , at which their stellar populations formed 
for (a) the critical model, Qo = 1, = and (b) a fiat model 
with Qo = 0.3, Ha = 0.7. In both plots, upper solid line cor- 
responds to 2f = 2, the upper dashed line to Zf = 3, the upper 
dot-dashed line to Z( = 5, the lower solid line to Z{ = 10, the lower 
dashed line to Zf = 20 and lower dot-dashed line to Z{ = 100. In 
both plots it is assumed that the stars are formed in a 1 Gyr 
star burst during which the star formation rate is constant. Solar 
metallicity is assumed. The mean values of the observed values 
of AK in equal cosmic time intervals of O.I/Tq are plotted as in 
Fig. ^, with error bars corresponding to the standard error in the 
mean of the K-magnitudes of both sets of observational data. 



4.2 Star Formation Epoch 

The epoch at which the stellar populations were formed was 
allowed to vary between redshifts of 2 and 100. The sense 
of the changes to the K-z relations are similar for all cosmo- 
logical models. Fig. ^ shows the K-band magnitude differ- 
ences of the models relative to the no-evolution track in the 
same cosmology for models with flo = 1, Qa = and f2o = 
0.3, Q.A = 0.7. The magnitude differences relative to the fidu- 
cial model are listed in Table 4 for Q,o = \, f^A = 0. Results 
have not been tabulated for models with f2o = 0.3, JIa = 0.7, 
as the magnitude differences within these models were not 
dissimilar to those in the Einstein-de Sitter cosmology. In 
Fig. P the values of AK are plotted against cosmic time as 
a fraction of To , where To is the present age of the Universe 
listed in Table 1 for Ho = 50kms"^Mpc"^ . 
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Table 3. Changes in the K magnitudes, AK, relative to the no-evolution model for four representative world models. The first part of 
this table gives the mean differences AK between the no-evolution model for these cosmologies and the observational data for all points 
within 20% of three different redshifts. The relations have been normalised to the mean observed magnitudes at small redshifts. Positive 
values indicate that the data points are brighter than the no-evolution track at that redshift. The second half of the table shows the same 
differences once account is taken of the passive evolution of the stellar content of the galaxies. For all four cases, the standard initial 
starburst is taken have the following parameters; Zf = 10, metallicity Z = Zq, Hq = 50kms~^Mpc~^, and the duration of the initial 
starburst 1 Gyr. Negative values indicate that the predicted evolution in K-band magnitude is less than that required by the data, that 
is, the modelled track is fainter than the mean of the observed relation at that redshift. 





Cosmological 




Sample 


AK from the no- 


-evolution track at different redshifts 




Model 






z=0.5 


z=1.0 


z=2.0 




= 1.0, Ua = 


0.0 


SCR Data 


0.597 ± 0.140 


1.022 ± 0.126 


1.780 ± 0.114 








6C Data 


0.706 ± 0.097 


0.433 ± 0.145 


1.136 ± 0.233 








Combined Data 


0.646 ± 0.089 


0.742 ± 0.115 


1.222 ± 0.210 


no 


= 0.1, Qa = 


0.0 


Combined Data 


0.778 ± 0.091 


1.110 ± 0.118 


1.945 ± 0.220 


no 


= 0.1, Oa = 


0.3 


Combined Data 


0.712 ± 0.091 


1.076 ± 0.119 


1.925 ± 0.220 


no 


= 0.3, Oa = 


0.7 


Combined Data 


0.789 ± 0.091 


1.101 ± 0.117 


1.758 ± 0.214 








Model 


AK difference between model and observations 


no 


= 1.0, Ca = 


0.0 


standard initial starburst 


-0.177 ± 0.089 


0.162 ± 0.115 


0.336 ± 0.210 


no 


= 0.1, Ca = 


0.0 


standard initial starburst 


-0.303 ± 0.091 


-0.262 ± 0.118 


-0.696 ± 0.220 


no 


= 0.1, ^A = 


0.3 


standard initial starburst 


-0.352 ± 0.091 


-0.310 ± 0.119 


-0.799 ± 0.220 


no 


= 0.3, SIa = 


0.7 


standard initial starburst 


-0.420 ± 0.091 


-0.295 ± 0.117 


-0.454 ± 0.214 



In Fig. |5|, it is assumed that the galaxies have solar 
metallicity. All models have been normalised to the same 
magnitude at the lowest redshifts, at which the mean K-z 
relation is well defined. The models diverge at earlier times. 
As expected, if the epoch of star formation is less than the 
reference value Zf = 10 the predicted tracks are more lumi- 
nous at redshifts z ~ 2 — 3. The largest differences occur in 
the cases in which the star formation is assumed to begin 
at Zf = 2 and 3, these being less than that at which radio 
galaxies are known to exist. Nonetheless, it is noteworthy 
that, unless the epoch of star formation were as unreason- 
ably low as Z{ = 2 — 3, changing the epoch of star forma- 
tion does not result in a particularly large change in AK. 
The reason for this can be understood from Fig. 3(a). By 
adopting 2.2 ^m as the wavelength at which the magnitude- 
redshift relation is determined, even at redshift 2 = 3, the 
changes in the luminosity of the galaxies at a rest wave- 
length of 500 nm are quite small and change monotonically 
with time. Indeed, the difference between observing the star- 
burst at ages 0.32 to 1.25 x 10^ years amounts to less than 
a factor of three at 5000A in the rest frame. This behaviour 
contrasts dramatically with what would be expected if the 
redshift-magnitude relation were determined in, say, the R 
waveband. Then, Fig. 3(a) shows that the variations in lu- 
minosity at rest frame wavelengths of 1500 A would be very 
sensitive to even small amounts of ongoing star formation. 
These calculations illustrate the considerable advantages of 
carrying out studies of the redshift magnitude relation in the 
K-waveband - the observations provide information more di- 
rectly about the bulk star formation activity of the galaxies. 

As expected, as Zf increases the predicted K-z relations 
become fainter at large redshifts and there is little difference 
between the tracks in both models if Zf is greater than 10. In 
the case of the model f2o = 0.3, f^A = 0.7, the fit to the data 
can be improved considerably if the star formation redshift 
~ 3 — 5 is adopted. Values of AK for the examples given 



in Fig. 5(a) at redshifts of 0.5, 1.0 and 2.0 are included in 
Table 4. 

4.3 The Value of Hubble's Constant 

Intuitively, it would seem that the choice of Hubble's con- 
stant would have a significant impact upon the form of the 
K-z relation since the cosmic time-scale is significantly dif- 
ferent if, say. Ho were are large as 100kms~^Mpc~^. Then, 
all cosmic time-scales would be reduced by a factor of two, 
allowing much less time for the evolution of the stellar pop- 
ulations of the galaxies and also observing them at very 
much earlier times at large redshifts. These changes, as well 
as those for Ho — 65kms~^Mpc~^, are given in Table 4. 
In the cases of models with larger values for Hq, an initial 
starburst of very short duration was adopted so that the 
galaxy is observed during the passive evolution phase; the 
tabulated values in Table 4 have included this change in star 
formation history. 

It is noteworthy that changing the value of Hubble's 
constant has a remarkably small effect upon the predicted 
evolutionary tracks. The reasons for this can be understood 
from inspection of Fig. 3(a). As was emphasised in Section 
3.2, so long as the analysis is restricted to redshifts less than 
3, the spectral evolution of the stellar energy distribution 
is monotonic with time and can be reasonably described 
by power-laws at different wavelengths from 0.5 to 2.2 fim. 
Thus, if the value of Hubble's constant is changed, the stellar 
mass of the galaxies required to normalise the K-z relation 
at low redshift is changed, but the differences in the values 
of AKov depend primarily upon the changing relative cos- 
mic times observed at different redshifts. In general, these 
are quite small and so changing the value of Ho does not 
greatly change the evolutionary corrections. If the redshift 
magnitude relation were determined at rest-frame UV wave- 
lengths, this would certainly not be the case. 
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Table 4. Changes in the K magnitudes, AK, relative to the standard passive evolution model in an Qq = 1, = world model for the 
representative world models. Results are given for redshifts z = 0.5, 1.0 and 2.0. Positive values of AK, which is defined as as Kq — Kmodeli 
indicate that the new model results in a greater increase in infrared luminosity as compared with the standard passive evolution model. 



Section 






AK from the observational data: 








z=0.5 


z=1.0 


z=2.0 




Standard Passive Evolution Model; 


-0.18 


0.16 


0.34 




o in o. nn -^^ in 7 7 

Si[) — l.U, — U.U, Zf — lU, Zj — ZjQ 










IGyr starburst, Uq = 50kms ^Mpc ^. 










Change from standard passive evolution model 


AK difference from standard passive evolution model: 






z=0.5 


z=1.0 


z=2.0 


(4.1) 


Oo = 0.1, Qa = 0.0 


-0.13 


-0.42 


-1.03 




Co = 0.1, Ua = 0.3 


-0.18 


-0.47 


-1.14 




Co = 0.3, f^A = 0.7 


-0.24 


-0.46 


-0.79 


(4.2) 


2f=2 


0.10 


0.27 


- 




2f=3 


-0.00 


0.10 


0.57 




2f=5 


-0.02 


0.04 


0.14 




2f = 10 


0.00 


0.00 


0.00 




2;f=20 


-0.04 


-0.05 


-0.14 




Zf=100 


0.00 


-0.02 


-0.12 


(4.3) 


J^O = 50kms-lMpc-l 


0.00 


0.00 


0.00 




Ho = 65kms~iMpc~l 


-0.02 


0.04 


-0.01 




Ho = lOOkms^^Mpc-i 


-0.01 


-0.05 


0.12 


(4.4) 


Z = 0.2Zq 


0.05 


0.22 


0.45 




Z = OAZq 


0.09 


0.10 


0.32 




Z = I.OZq 


0.00 


0.00 


0.00 




Z = 2.5Zq 


0.13 


-0.04 


-0.16 


(4.5) 


Instantaneous SFR 


-0.04 


-0.09 


-0.21 




IGyr Constant rate 


0.00 


0.00 


0.00 




Exponential SFR 


0.05 


0.07 


0.05 


(4.5) 


Extra starburst, 2.5% mass, lO^years old 


0.04 


0.07 


0.17 




Extra starburst, 2.5% mass, lO^years old 


0.32 


0.28 


0.22 




Extra starburst, 2.5% mass, lO^years old 


0.04 


0.04 


0.05 


(4.5) 


Extra starburst, lO'' years old, 0.5% mass 


0.07 


0.06 


0.05 




Extra starburst, lO'' years old, 1.0% mass 


0.14 


0.12 


0.09 




Extra starburst, lO^years old, 2.5% mass 


0.32 


0.28 


0.22 




Extra starburst, lO^years old, 5.0% mass 


0.56 


0.50 


0.40 



4.4 Varying the Metallicity of the Galaxies 



The metallicity of the galaxies can significantly influence 
their evolutionary tracks. Rather than attempting to adopt 
some arbitrary chemical history, the evolutionary codes have 
been run for metallicities of 0.2^©, OAZq, Zq and 2.5Zq. 
The results are displayed in Fig. ^ for the same two cos- 
mologies discussed in Section 4.2 and for Z{ = 10. For the 
Einstein-de Sitter cosmology, these results are also tabu- 
lated in Table 4. All four K-z tracks are normalised to the 
same magnitude scale as the fiducial model, rather than to 
the same magnitude at the present epoch. This procedure 
has been adopted so that at all times the variations due to 
changes in metallicity are clearly displayed. The sense of the 
changes is that, the lower the metallicity, the more luminous 
the galaxy at early epochs. It is apparent from Fig. ^(b), that 
assuming a lower metallicity at large redshifts improves the 



fit of the Slo ~ 0.3, SIa = 0.7 model as compared with the 
fiducial model. 

The lowest metallicity track is initially the brightest 
(largest AK), and the 2.5^© track the faintest, but the 
tracks cross over at low redshift. The greatest variation be- 
tween the tracks occurs at early epochs and, as we have seen 
from Section 3.2, it is at these young ages that the galaxies 
undergo the greatest amount of evolution. The early evolu- 
tion of the spectral energy distribution is highly sensitive to 
the metallicity of the galaxy because, for a star of a given 
mass, changing the metallicity leads to a change in the age 
at which it moves from the main sequence onto the AGB; 
for lower values of Z/Zq the stars evolve more rapidly since 
there is less line blanketing and lower opacity. This has two 
consequences: firstly, the more rapid evolution means that 
lower metallicity stars were more luminous at early times; 
secondly, the stars were also shorter lived and at later times 
there are fewer stars present in the older galaxies, which 



© 2001 RAS, MNRAS 000, §-0 



K-z Relations for Luminous Radio Galaxies 11 




1 1 1 1 1 — 

0.5 



Cosmic Time as a fraction of To 





V ■ ■ 




(N - 










▲ 


O - 







1 1 1 1 1 — 
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Cosmic Time as a fraction of To 
Figure 6. The values of AK as a function of cosmic time relative 
to the no-evolution model for different assumed metallicities for 
(a) the critical model, = l,f^A = and (b) a flat model 
with Ho = 0.3, Ca = 0.7. The AK-cosmic time relations for the 
different metallicities are as follows: dotted line Z = 0.2Zq \ dot- 
dashed line Z = OAZq; solid line Z = Zq and dashed line Z = 
2.5.^0. The models involve a 1-Gyr constant starburst, in which 
star formation begins at Z{ = 10. Symbols are as in Fig. ^. 

therefore appear fainter than a higher metallicity galaxy of 
a similar age. 

It is interesting that the divergence of the different 
metallicity tracks occurs when the Universe was roughly half 
its present age, corresponding to the redshift at which the 
K-z relations for the SCR and 6C samples begin to diverge. 
It is just conceivable that such an effect could account for the 
difference between these samples, but it would require those 
6C galaxies which are less luminous than the SCR galax- 
ies to have higher metallicities. There is no obvious reason 
why greater radio luminosity should be correlated with lower 
metallicity of the stellar populations of the radio galaxies. 

The diagrams in Fig. M can be used to construct ap- 
proximate evolving galaxy models in which the metallicities 
of the galaxies increase with cosmic epoch. Such evolution 
tracks can be found by changing from one locus to another 
as the metallicity increases. The range of models shown in 
Fig. ^ illustrates how such changes would alter the predicted 
K-z relations. 



4.5 Star Formation History 

The results so far have involved an initial starburst during 
which the star formation rate has been assumed to be con- 
stant. It is important to consider also models in which there 
is ongoing star formation, and in particular those in which 
the occurrence of a radio source event might be associated 
with a major star formation event in the host galaxy. This 
could occur if, for example, the central black hole were fu- 
elled following some form of galaxy merger. To understand 
the significance of these results, we need to look more closely 
at the origin of the K-band flux density. 

The K-band luminosity of evolving galaxies arises from 
two stellar populations. Stars with mass M>2M0 evolve be- 
yond core and shell helium burning to become asymptotic 
giant stars; those with M<2M0 become red giant stars form- 
ing the old red giant branch seen, for example, in globular 
clusters. For the present study, we are mainly interested in 
the evolution of the overall stellar population of the galax- 
ies over cosmic time-scales greater than lO^years. This age 
corresponds to the main sequence lifetime of a star of mass 
2M0. Therefore, at ages greater than lO^years, the evolution 
is dominated by stars evolving onto the standard red giant 
branch and the arguments given in section S.2 can account 
for the decreasing luminosity of the galaxy. On the other 
hand, if there are epochs of rapid star formation, the bulk 
of the near-infrared luminosity is associated with stars on 
the asymptotic giant branch. For the most luminous stars, 
evolution onto the asymptotic giant branch takes place af- 
ter a few times lO^years. Figure 0(b) shows dramatically the 
large increase in K-luminosity of an instantaneous starburst 
when it is about lO^years old associated with the formation 
of AGB stars. On the other hand, if the starburst is of longer 
duration, say lO^years as illustrated in figure |^(c), the effect 
is somewhat diluted by the continued formation of all types 
of massive stars. 

Considering first variations in the long-term star forma- 
tion history of the galaxy populations. Fig. |^ displays three 
different types of star formation history: (i) the standard 1- 
Gyr starburst, (ii) an instantaneous starburst of negligibly 
short duration, and (iii) a model in which the star forma- 
tion rate decays exponentially with an e-folding time of 10^ 
years. There is very little difference between the three mod- 
els, all three of which have been normalised to the same 
stellar mass at late cosmic times. 

The brightest model (largest AK) except at the very 
earliest times is the exponential star formation model, and 
the faintest the instantaneous burst model. These results are 
not surprising since, it is only the behaviour after lO^years 
which is relevant to the K-z relation for z < 3. Once star 
formation ceases, the continued star formation of the ex- 
ponential model, and sustained production of AGB stars 
(albeit in small numbers), is sufficient to make this model 
the brightest track at low redshift. The most evolved model, 
the instantaneous starburst, has the lowest luminosity at 
low redshift. These differences between the models are very 
small, as can be seen from the values given in Table 4. 

Fig. ^ shows the standard model coupled with the SEDs 
of instantaneous starbursts at three different ages (10®, lO'' 
and 10* years old), the additional starburst accounting for 
2.5% of the total stellar mass of the galaxy. The model with 
the additional burst occurring the longest time before ob- 
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Cosmic Time as a fraction of To 
Figure 7. The values of AK as a function of cosmic time relative 
to the no-evolution model for the critical model, Qq = 1, Oa = 
0, for galaxies modelled with different star formation histories. 
Dashed line; instantaneous starburst. Solid line: 1 Gyr starburst 
with constant star formation rate. Dot-Dashed line: Exponentially 
decaying starburst with 10® year e-folding time-scale. Galaxies are 
modelled assuming solar metallicity, and Zf = 10. Symbols are as 
in Fig. |{a). 

servations, 10* years, gives the least change from the pre- 
dictions of the underlying model. Similarly, the 10^ year old 
recent burst also results in only a slight difference at low red- 
shift. In the first case, the additional young (10* years old) 
high mass stars have evolved beyond the AGB, and in the 
second case (lO'' years old) they have not yet had enough 
time for their evolution to reach that stage. Fig. ^ shows 
that the rest frame K-band fiux of an evolved instantaneous 
burst at these ages is much closer to that for the underlying 
star formation model than the 10^ year old track. At high 
redshift, however, we need to consider a somewhat bluer 
rest-frame wavelength region. Whilst there is a small dif- 
ference for the 10* year old model, the 10^ year old recent 
burst is slightly brighter still, which is seen in Fig. ^. The 
model with the lO'^ year old recent burst component shows 
the greatest change from the standard model, as expected 
from Fig. ^ The AK values for these models given in Table 4 
are solely due to the additional recent burst component. The 
additional burst consistently brightens the K-z tracks at late 
times. The reason for the effect being most pronounced at 
small redshifts is that the luminosity of the starburst is as- 
sumed to be independent of redshift. Since the K-luminosity 
increases in roughly power-law fashion with increasing red- 
shift, the luminosity of the starburst is soon swamped by 
the majority population of red giant branch stars. Fig. 
displays the results for the same set of models as in Fig. rT 
but the mass fraction of the additional starburst is varied 
rather than its age. The tracks in this figure involve an ad- 
ditional 10^ year starburst and four mass fractions. 

If the mass of the starburst as a function of redshift 
remains constant, the overall shape of the K-z relation is 
flattened at late times and does not result in an improved 
flt to the observational data. Some scatter may be intro- 
duced into the K-z relation, because the radio galaxies are 
observed at a range of ages from lO'' to 10* years. In or- 
der for starbursts to change significantly the shapes of the 
K-z relations, the additional starburst mass would need to 
change with redshift. This is quite possible, as there is more 
gas available in any merger- induced starburst at high red- 
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Cosmic Time as a fraction of To 
Figure 8. Changing AK with age of the universe, for galaxies 
modelled with complex star formation histories with an addi- 
tional starburst occurring a fixed time before the observations, 
as outlined in section 2.2. The plotted tracks represent the fol- 
lowing models: solid line - no additional starburst, dashed line - 
additional starburst occurring 10^ years prior to the observations, 
dot-dashed line - starburst 10^ years old, dotted line - starburst 
10* years old. Galaxies are modelled with solar metallicity, and 
Z[ = 10, in an Einstein-de Sitter universe. The extra starburst 
accounts for 2.5% of the total stellar mass of the galaxy. Symbols 
are as in Fig. ^. 

shift. The same would be true if the starburst were radio jet 
induced (e.g. Rees 1989) , which provides a possible explana- 
tion of the alignment effect (McCarthy et al. 1987, Chambers 
et al. 1987). Figs. ^ and 9 can be used to estimate how sig- 
nificantly the K-band luminosity could have changed if the 
magnitude of the starburst were to vary with cosmic epoch. 
Just as in the case of evolving metallicity, so in the case 
of evolving starbursts, the evolutionary history can be esti- 
mated by moving from one track to another in Figs. |^ and 
9. For example, in Fig. ^, if the early starbursts involved 5% 
of the galaxy mass, the evolutionary track would begin on 
the dotted line; if this subsequently decreased to 2.5% and 
1%, the track would change to the dot-dashed and dashed 
lines respectively. A 5% starburst by mass is the absolute 
maximum allowable given the constraints of the observed 
ultraviolet luminosities for the SCR radio galaxies at z ~ 1 
(e.g. Best et al. 1998). In this extreme case the magnitude 
increase in K is at most about 0.5 magnitudes. 



5 DISCUSSION AND CONCLUSIONS 

The above analyses of the K-z relations for a wide range of 
plausible cosmological models have established a number of 
useful general rules concerning the potential use of the lumi- 
nosities of massive galaxies as standard candles. The value 
of carrying out these studies in the K- waveband at 2.2 
for redshifts z < 3 is clearly illustrated by the synthesised 
spectral energy distributions shown in Fig. ^(a) in which it 
can be seen that the evolution of the dominant stellar popu- 
lations can be predicted with some confidence. As has been 
discussed by previous authors, the dominant evolutionary 
trends are associated with the evolution of stars in the mass 
range 1-2 Mq from the main sequence onto the red giant 
branch. An equally important result is that, over this range 
of redshifts, the choice of cosmological model is much more 
important than the corrections associated the evolution of 
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Cosmic Time as a fraction of To 
Figure 9. Changing AK witli age of tlie universe, for galaxies 
modelled with star formation histories similar to those in Fig. ^ 
The plotted tracks represent models with the following stellar 
mass fractions involved in the extra starburst: solid line - 0.5%, 
dashed line - 1.0%, dot-dashed line - 2.5%, dotted line - 5.0%. 
Galaxies are modelled with solar metallicity, and Zf = 10, in an 
Einstein— de Sitter universe. The extra starburst is modelled as 
occurring lO'^ vears before the observations were made. Symbols 
are as in Fig. H. 

the stellar populations of the galaxies, as is vividly illus- 
trated by the values in Table 2. Thus, provided standard 
galaxies can be identified over the redshift range < z < 3, 
there are good prospects for obtaining useful information 
about the choice of cosmological model. 

The problem with this programme is the identification 
of such a class of standard galaxy, particularly in the hght 
of the favoured picture of hierarchical clustering to account 
for the formation of structures such as galaxies and clus- 
ters of galaxies over the same redshift interval. The analysis 
of Best et al. (1998) showed that the SCR radio galaxies 
seemed to form a population of galaxies of constant mass 
within the framework of the fio = l,fiA = world model, 
although they described a number of reasons why this was 
a somewhat surprising result. The present calculations have 
extended that analysis to a much wider range of cosmolog- 
ical models and considered in more detail the effects of the 
evolution of the stellar populations of the galaxies. The re- 
sults shown in Fig. ^ and documented in Table 2 show that 
all the new models are a significantly poorer fit to the large 
redshift data. In particular, the principal cause of the poorer 
fits is the choice of cosmological model. Expressed in phys- 
ical terms, in the critical f2o = l,f^A = model, galaxies 
at a given redshift are observed at significantly smaller dis- 
tance measures than in all the other models considered and 
so have the brightest K-magnitudes. In all the other mod- 
els considered, a literal interpretation of the results would 
be that the radio galaxies at redshifts z > 1 have greater 
stellar masses than the radio galaxies nearby. The analysis 
has shown that this result is independent of the choice of 
Hubble's constant. 

It is of particular interest to find ways in which the 
favoured flat cosmological model with Q.o = 0.3, = 0.7 
can be reconciled with the observational data. As the results 
in Table 4 indicate, to produce an improved fit to the ob- 
servational data one or more of the following changes to the 
standard evolution picture must be introduced: 

• A lower redshift of star formation, 3 < < 4, pro- 



vides a much improved fit to the observed K-z relation for 
this cosmological model (Fig. ^). A problem with this im- 
proved fit is that the galaxies would then be expected to 
be very much more luminous in optical wavebands with sig- 
nificant evidence for younger stellar populations, whilst the 
red colours of distant radio galaxies suggest a much larger 
luminosity-weighted stellar age (e.g. cf Dunlop et al., 1996). 
Further, the observation of well-formed galaxies which are 
powerful radio sources out to redshift z = 5.19 (van Breugel 
et al., 1999) indicates that at least some of these galaxies 
formed at higher redshifts. 

• The metallicities of the galaxies may have been signif- 
icantly less at redshifts z ~ 2 than they are at the present 
epoch. As shown in Fig. 6(b), decreasing the metallicity of 
the galaxy to only about 20% of the present cosmic value 
improves significantly the fit of the model to the observa- 
tions. 

• Starbursts associated with the radio source events can 
lead to a significant enhancement of the K-band luminosi- 
ties of the galaxies, due to the rapid evolution of luminous 
massive stars onto the asymptotic giant branch. This could 
potentially improve the fit of the models in lambda cosmolo- 
gies if associated starbursts brighten the distant galaxies, 
but not those at redshifts less than 1. This is quite plausible 
since such starbursts could be enhanced at high redshifts 
due to the extra gas available at these epochs for forming 
stars, either in a galaxy merger or triggered by the radio jet; 
indeed additional blue luminoisity is observed in powerful ra- 
dio galaxies at high redshifts (the alignment effect). There 
is clearly much scope for arbitrarily adjusting the masses, 
timescales and cosmic evolution of the starbursts to improve 
the fit to the K-band data. However, the excess blue lumi- 
nosity associated with the alignment effect requires only a 
few xIO^Mq of newly formed stars (< 5% of the galaxy 
mass), even were it all to be associated with star formation 
activity (e.g. Dunlop et al. 1989; Best et al, 1998): this sets 
a strong constraint upon the luminosities and durations of 
the bursts. 



The alternative approach is to accept that, even after 
accounting for stellar evolution, the radio galaxies at z > 1 
may indeed be more luminous, that is, have a greater stellar 
mass, than those at small redshifts. This conclusion would 
mean that the low and high redshift radio galaxies would not 
form a single evolving population, and would be surprising in 
terms of hierarchical clustering scenarios according to which 
such massive galaxies at high redshift are extremely rare. 
Given the growing evidence that the distant radio galaxies 
belong to cluster environments, it may be plausible that they 
had undergone more mergers than the nearby radio galaxies, 
which are found in more rarified environments: this possibil- 
ity could be tested by carrying out cross-correlation analyses 
of radio galaxies of different K-luminosities with field galax- 
ies, although this becomes a demanding project for the ra- 
dio galaxies at z > 1. The similar characteristic sizes of the 
SCR radio galaxies at redshifts 1 and 0, however, sets con- 
straints upon the intrinsic differences in luminosities which 
effectively rule out the most extreme cosmological models, 
those with To x ffo ^ 1- 
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